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Abstract 

It is shown that the interaction of large amplitude electromagnetic waves with a hot electron- 
positron (e-p) plasma (a principal constituent of the universe in the MeV epoch) leads to a bunching 
of mass, energy, and angular momentum in stable, long-lived structures. Electromagnetism in the 
MeV epoch, then, could provide a possible route for seeding the observed large-scale structure of 
the universe. 
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I. INTRODUCTION 



It is widely believed that the currently observed large-scale structure of the universe 
(the clusters and superclusters of galaxies) grew gravitationally out of small density fluc- 
tuations ffl. The imprint of this density variation in early universe was left on the cosmic 
microwave background radiation in the form of spatial temperature fluctuations. The grav- 
itational origin, however, can be only a part of the story; the gravity can enhance, but it 
can not produce these fluctuations. The quest for the physical process(es) which produced 
the initial matter-density fluctuations, has led to the emergence of the following two leading 
mechanisms: inflation, and topological defects P, 0]. Of these the former is, perhaps, the 
most thoroughly investigated. According to this mechanism there exists, in the evolution 
of the universe, an early inflationary period in which the universe expands so rapidly (ex- 
ponentially) that quantum fluctuations get trapped in the expansion. By the end of the 
inflation, therefore, small irregularities covering a wide range of length-scales permeate the 
entire universe. Gravitational instability then acts on these small initial irregularities and 
enhances the concentration of matter from which galaxies and clusters of galaxies eventually 
emerge |3J. Theory of such processes is, by no means, complete and much needs to be done 
to determine whether tiny quantum fluctuations can provide a strong enough template for 
gravitational condensation to finally create the structures that we observe today. Cosmolo- 
gist, usually, rely on the speculated existence of non-baryonic dark matter to augment the 
gravitational force to aid and accelerate the structure formation. The question is far from 
settled! 

It is natural, then, to look elsewhere for the source of the seed density "fluctuations." 
An obvious possibility is to explore if electromagnetic interactions taking place in a plasma 
(known to be the source of a whole variety of linear as well as nonlinear waves) can cause 
the required density perturbations. In the standard cosmological model of the hot Universe 
(the Big Bang model), it is estimated that temperatures as high as T ~ 10 10 K ~ 1 MeV 
prevail up to times ~ 1 second (t ~ 1 sec) after the Big Bang. In this epoch, the main 
constituents of the Universe are photons, neutrinos and anti-neutrinos, and e-p pairs ^ |J. 
As the plasma cools down the annihilation process e + + e~ — » 7 + 7 dominates and the e + e~ 
pair concentration goes down. Since the equilibration rates are fast in comparison with the 
changes in plasma parameters, an equilibrium e-p plasma should be present in the MeV 
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epoch of the early Universe. It is this plasma-dominated era in which we will seek the seeds 
for future structure formation. 

Relativistic e-p plasmas have been investigated quite extensively ||. Tajima and 
Taniuti || suggested that collective processes in these plasmas could lead to interesting 
consequences for structure formation. In the e-p plasma of the early Universe, localized low 
frequency electromagnetic (EM) waves (%u <C T) could propagate as an envelope soliton 
due to the interaction with sound waves. Plasma density variations related to these solitons 
could potentially be useful towards structure formation in the Universe. However the anal- 
ysis H is based on a one dimensional formulation, and the corresponding soliton solutions 
are likely to be unstable in higher dimensions. Berezhiani and Mahajan |7| argued that in 
the MeV epoch of the Universe, although the e-p pairs form the dominant constituent of the 
plasma, a minority population of heavy ions is also present due to the baryon asymmetry. 
They were able to show that, under appropriate conditions (when the plasma is transparent, 
i.e., uj u e , where uj (u> e ) is the pulse (plasma) frequency), the resulting e-p-ion plasma 
supports the propagation of stable, nondiffr acting and nondispersing EM pulses (light bul- 
lets) with large density bunching. It was further shown in Ref. [|J that these bullets are 
exceptionally robust: they can emerge from a large variety of initial field distributions, and 
are remarkably stable. Note that the characteristic dimensions of such matter-filled light 
pulses are proportional to the electron to ion density ratio, and tend to be considerably 
larger than the skin depth (A = c/uj e ). The implication is that when one deals with the EM 
structures whose characteristic dimensions (of the spatio-temporal inhomogeneities) are of 
the order of the skin depth, the baryon asymmetry affects can be safely neglected, and the 
dynamical system can be assumed as a pure electron-positron plasma. 

In the present paper we examine the propagation of strong EM radiation in a hot pure 
e-p plasma with the explicit aim of finding soliton-type solutions. The plasma is assumed 
to be transparent. We demonstrate that the dynamics of the EM field envelope is governed 
by a generalized nonlinear Schrodinger equation (NSE) with a defocusing nonlinearity. In 
ID, this equation admits dark soliton solutions, while in 2D, the so called vortex soliton 
solutions are also possible. 

Dark solitons exist as dips on a continuous-wave background field. Being stable in one 
dimensional, they appear as dark-stripe solitary waves in bulk medium. However, such 
stripes are unstable to transverse modulations that results in the induced generation of 
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vortices with alternating polarities. 

Vortex solitons, which are the most fundamental 2D soliton solutions of NSE with an 
angular 2tt phase ramp, appear as local dark minima in an otherwise bright background. 
Vortex solitons have been recently observed in materials with a defocusing optical nonlinear- 
ity — the dynamics of laser beams in these materials is generally described by the NSE || . 
Since the electromagnetic vortices (EMV) carry angular momentum that is conserved dur- 
ing propagation, the generation of Vortex solitons in an e-p plasma is a potent mechanism 
for creating domains with definite angular momenta even out of an initial field distribution 
devoid of angular momentum. To keep the total angular momentum zero, domains of equal 
and opposite angular momentum must be created in pairs. 



II. FORMULATION 



We use the following set of relativistic hydrodynamic equations in dimensionless form 



+ ^VP ± = T^± [v* x (V x A)] =f V0, (2) 
^ + V ■ (n^) = 0, (3) 

along with the field equation (in the Coulomb gauge V • A = 0) 

cP 1 

—A-V 2 A+—V(j)+(n-v--n + v + ) = 0, (4) 
at 1 at 

where p^ 1 = 7 ± t> ± with factor 7 =b = [1 + (p ± ) 2 ] 1/ ' 2 , d±/dt = d/dt + v± ■ V is the comoving 
derivative, and G ± = K^(l/T ± )/K2(l/T ± ) with K n the n-th order modified Bessel functions 
of the second kind. The superscript labels the particles, electrons (— ) and positrons (+) 
respectively. In these equations the time and space variables are in units of the electron 
plasma frequency to e = (47re 2 n /m e ) 1//2 and the collisionless skin depth c/uj e respectively, 
the field potentials (0, A) are in the units of m e c 2 /e, and the relativistic momentum vector 
p ± is in units of m e c. The particle number density n is normalized by the equilibrium 
density n = and plasma temperature (T^ is measured in units of m e c 2 . The pressure 
= n^T 1 * 1 , where nf is the density in the rest frame of the fluid element (nf = n ± /7 ± ). 
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The function G(z) defines the "effective" temperature dependent mass of the particles and 
has the following limiting expressions: G ~ 1 + 5/2z for z ^> 1, and G ~ 4/z if z <C 1. 
From Eqs. ([j])-(|3]) it is straightforward to derive the adiabatic equation of state flT0|| : 



n±/T 



± 



exp(-G ± /T ± ) = const., (5) 



7 ±K 2 (1/T±) 

which, at nonrelativistic temperature (T ± <C 1), reduces to the standard adiabatic relation 
(nf / {T ± )' i l 2 = const) for an ordinary gas. In the ultrarelativistic limit (T 1 * 1 ^> 1), as 
expected, Eq. (||) describes the photon gas (n^/fT 1 * 1 ) 3 = const). In the ultrarelativistic 
case, one should take into account the radiative pressure Pr = aT 4 (a = n/45h 3 c 3 ). For 
simplicity we neglect this less important effect for the current considerations. Notice that 
in the MeV epoch, the plasma temperature m m e c 2 (i.e. z « 1), and G ~ 4 leading to 
an effective mass of e-p pairs to be m c s ~ 4m e . Since the particle masses are just a few 
times larger than their rest mass at these temperatures, the e-p plasma can be considered 
as a two component fluid rather than a photon gas. 

We consider the propagation of circularly polarized EM wave with a mean frequency 
u, and a mean wave number k along the z axis. The choice of circular polarization is not 
restrictive, it simplifies the analysis by preventing harmonic generation. The vector potential 
can be represented as: 

A± = —(x + \y) A±(r±, z, t) exp(i&;,2 — iut) + c.c, (6) 

where A± is a slowly varying function of r and t (k ^> V, u ^> d t ). The unit vectors 
x and y define two mutually perpendicular axes in the plane normal to the direction of 
wave propagation. The Coulomb gauge condition leads to the relation A z = (i/k)(V± ■ 
A±) A±. Consequently the effects related to A z will turn out to be negligibly small. We 
shall now follow standard methods to analyze the system. In the slowly varying amplitude 
approximation, the transverse, high-frequency component of the equation of motion yields 
the simple relation between the particle momentum and the vector potential 0, 

p±G ± = T A ± . (7) 

The low frequency motion of the plasma is driven by the ponderomotive pressure 
( ~ (p^ 2 ) °f the high-frequency EM field, and it does not depend on the sign of the 
particles' charge. If we assume that in equilibrium the electron and the positron fluids 
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have equal temperature (T^ = T ), their effective masses will also be equal (G^ = G), 
and the radiation pressure will impart equal low-frequency momenta to both fluids allowing 
the possibility of overall density changes without producing charge separation. The charge 
neutrality conditions n~ = n + = N, <fi = will be assumed in the rest of the paper. It 
is also evident that the symmetry between the two fluids keeps their temperatures always 
equal (T = T) if they were equal initially. 



Considerable simplification results when we invoke the wide beam approximation |TT 
we assume that the longitudinal variation of the field envelope is much stronger than the 
transverse variation, i.e. L z , the characteristic length along the propagation direction, is 
much shorter than L±, the characteristic length in the transverse plane. This approximation, 
coupled with charge neutrality, allows us to extract from Eqs. (|l])-(@)> the following, leading 
order description for the low frequency response: the equation of motion 

d l^NT = _J_d^l 

dt P Ndz 7 " 2 7 G dz ' { ) 

and the "energy" conservation equation: 

i^-i^ = JL?ky! (9 ) 

dt 1 Ndt 7 2 7 G dt ' { } 

Here we have used the condition that the ponderomotive pressure gives equal longitudinal 
momenta to both electrons and positrons (pf = p). Notice that the assumed circular 
polarization of the EM field insures that the relativistic factor 7 does not depend on the 
"fast" time (1/cu) scale; it can be written as 
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1 + L -^t+p 2 



1/2 

(10) 



Substituting Eqs. (^) and (0) into Eq. (||), we find that the slowly varying amplitude A± 
must satisfy 

2icu(d t + v g d z )A ± + VlA ± 

2N 

+ (d 2 z - d?)A± + ^ 2 ~ k 2 )A± ~ -qA ± = 0, (11) 

where v g denotes the group velocity of the carrier waves, v g = k/u. 

We are still not quite done with simplifying assumptions. We seek solutions which vary 
slowly with time in a frame comoving with the wave, that is, in the frame propagating with 
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the group velocity v g . The transformation £ = z — v g t, t = t with the condition v g d^ 3> d T 
helps implement this approximation. Equations (|8|) and (0) can now be combined to derive 

§^[G(l - v g p)] = 0, (12) 

the implied constant of motion is to be determined from the boundary conditions. We 
demand p and A± to be zero at infinite £, but allow them to be finite as r± — > oo. Integrating 
Eq. ( |T2"D leads to (T is the particle temperature at infinity) 

G{T){~f-v g p)=G (T ), (13) 

which is readily solved for an explicit expression for the longitudinal momentum in terms of 
the transverse vector potential 



2 G 

P = v glg — 



1 -^ (7 ' G °- G2 -'^ |2)1/2 



(14) 



where 7 g = 1/(1 — v 2 ) 1 ^ 2 is the "effective relativistic factor" associated with the group 
velocity of the wave; it is not to be confused with the particle 7. The continuity equation 
can be similarly integrated to determine the particle density (after using Eq. (|14]) for p): 

7 G {^Gl-G 2 -\A ± \ 2 y/ 2 ' y } 

Substituting Eq. ([TJ) into Eq. (|TTD , we obtain the following nonlinear Schrodinger equa- 
tion for the complex amplitude A±, 

2iud T A ± + VlA± + -jdlA± 
^9 

where the wave frequency u satisfies the dispersion relation uj 2 = k 2 +2/Go, implying that the 



parameter 7^ = uyG /2 (7^ = (u/uj e )^G /2 in physical quantities). The set of Eq. 
and the equation of state (|5|) (in which the relation (|15|) could be easily incorporated) 
constitutes a complete description of the dynamics of strong EM waves in relativistic e-p 
plasma in the wide beam approximation. 



We remind the reader that Eq. (|T6|) was derived under the assumption ^> Vj_ (i.e. L z <C 
L±). In spite of that, for a highly transparent plasma (j g ^> 1), the second, "diffractive" 
term can be the same order or even greater than the third, "dispersive" term. For this paper, 



we will not attempt the general solutions of this quite complicated set of equations; we will 
simply deal with waves for which the plasma is so highly transparent that the diffractive 
term dominates. Using 7 9 G 3> G, and neglecting the dispersive term, the NSE simplifies 
to 



id T A± + -V?A ± - 2 



12 \ -1/2 



A ± = 0, (17) 



where the following renormalizations are used: t/2luGq — ► r, r±/\ / 2Go — > rj_. 

The vector potential \A±\ is restricted from above by the condition \A±\ < 7 s G*o- This 
restriction is necessary for the validity of the hydrodynamic treatment for the particles. For 
larger amplitudes, the electromagnetic waves are overturned causing multistream motion 
of the plasma requiring a kinetic description. Notice that despite the upper bound on the 
amplitude of the vector potential, the EM field can be still relativistically strong, i.e., the 
normalized \A±\ 3> 1, since j g 1. 

III. STATIONARY SOLUTIONS 

In the NSE derived above the diffractive and nonlinear terms have opposite signs and as a 
consequence Eq. fllTD does not admit transversely localized solutions (also called the bright 
solitons). Any localized initial EM field, therefore, will undergo transverse spreading during 
propagation. The NSE with a defocusing nonlinearity can, however, support stationary 
structures with asymptotically (at infinity) nonvanishing fields. Dark solitons in ID, and 
vortex solitons in 2D, are the fundamental representatives of such solutions. In the extreme 
low amplitude limit, \A±\ <C 7 9 Go, Eq. (|I7]) reduces to the NSE with a cubic nonlinearity. 
In one- dimensional geometry we have: 

dA± ld 2 A ± 1 . , , 2 A 

^ + 2^-^Gl^ A - = °' (18) 



This equation is exactly integrable via the inverse scattering method [12]|, and its one-soliton 
solution can be written as HI 



A ± (x, r) = 7 5 G ^o(« tanh 6 + i(3) e~ iA ° T , (19) 
where 

6 = aA (x -0Aqt). (20) 



S 



Here Aq is a measure of the asymptotic fields at the spatial infinity and a and (3 are con- 
stants with a 2 + f3 2 = 1. The solutions with a nonzero value at the center of the dip, is 
termed the 'gray soliton' to distinguish it from the 'black soliton' (zero amplitude at the dip) 
corresponding to /3 — 0. The dark solitons of this class of NSE's do not have any threshold 
values for their excitation unlike bright solitons (of the appropriate equations) which do. In 
other words, dark solitons can be created by an arbitrary small initial dip on a homogeneous 
background. 

In two transverse dimensions (2D), a dark soliton represents a dark stripe imposed on a 
homogeneous bright background. It is well-known that such a stripe is unstable to transverse, 



long wave length modulations |T3| . The instability causes the stripe to split into a sequence of 
Vortex solitons of alternative polarities. The vortices are dark holes on a bright background, 
with a nested phase dislocation of the order m = ±1, ±2, ... at their core. 



The Vortex soliton solutions of NSE were first suggested by Pitaevskii |L4j as topological 
excitations in an imperfect Bose gas in the superfluids. The ability of some electromagnetic 
systems (like the e-p plasma) to simulate fluid dynamical phenomena (like vortex formation) 
can be demonstrated by applying the Madelung transformation: A± = ^fp exp(h/>) to the 
defining equations [fOJ. The transformation converts the original set to one that is similar 
to the fluid hydrodynamic equations with a fluid "density" p, and fluid "velocity" v = V±i/j. 
Vortices can exist despite the potential nature of the "fluid" flow. Indeed, the Madelung 
transform is singular at the points where p = 0; these are just branch points where the real 
and the imaginary parts of the field become zero, while the velocity circulation § v ■ dl = 
27rm, where the integration is done on a closed path enclosing the singular point, and the 
integer m is known as the topological "charge" of vortex. Thus the vortex soliton is a 
topological structure; it can disappear only when annihilated by a vortex soliton of the 
opposite "charge." The development of the transverse instability of a dark soliton has close 
parallels in hydrodynamics: for instance the Kelvin-Helmholtz instability, which occurs 
when the boundary between two flows develops the so-called vortex streets ||15|| . Since a 
dark solitary stripe does not carry any topological "charge," it is evident that vortices have 
to be born with equal and opposite topological "charges." 

It is straightforward to show that e-p plasmas can support large amplitude dark solitons 
as well. In the general case (amplitude large, but subject to the condition \A±\ < 7 5 Go), 
we cannot construct analytic solutions even in ID. It is possible, however, to extract the 
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general properties of the solution by using reasonably simple techniques especially when A± 
has the time dependence 

A(x, t) = A(x) exp(-iAr), (21) 

where A = A^/'jgGo is the normalized amplitude, and is always less than unity. Here A is 
so-called nonlinear frequency shift. This time-dependence implies that the amplitude square 
is stationary (what follows, therefore, are classed as stationary solutions), and the dip of the 
wave does not propagate in the comoving frame with quite the group velocity of the linear 
wave. 

The ID Eq. (|18|) now can be cast in the form 

^A + V '(A) = 0, (22) 
where the prime on V denotes the derivative with respect to A, and 



V{A) = {X + 2)A 2 + A\/1- A 2 -A, (23) 



denotes the potential. The resemblance of Eq. ( |22|) to the one obeyed by a Newtonian 
particle in a nonlinear potential suggests an obvious method for analysis. One can easily 
prove that bounded solution exists provided nonlinear frequency shift is positive (A > 0). 
The profile of the potential, shown in Fig. |T| for A = 1, reveals that the dark soliton solution 
may reside in the potential well. Equating d 2 A/dx 2 to zero, we estimate the upper bound 
on A 



2 



A '^f-{—2)- (24) 
The lowermost value of the amplitude \A\ is zero, that is, we recover the dark soliton. Note 

that for small values of A ( <C 1), A u ^ — > 0, while A u ^ — > 1 for A 1. 

Similarly one can show the existence of vortex soliton solutions in 2D. We shall again 

seek stationary solutions in 2D polar coordinates (r, 9). The ansatz 

A = A(r) exp(im# - iAr), (25) 



with A(r) real, and with the perpendicular Laplacian operator given by 

d 2 Id m 2 



converts Eq. (|17D to the ordinary differential equation (similar to Eq. (p2|) ), 

^ 2 A v'(A)- m% A (27) 

dr 2 r dr r 2 ' 

10 




FIG. 1: The form of the potential V{A) in Eq. (ph is shown for A = 1 



where the potential V(A) is the same as the one- dimensional expression given by Eq. (f23"|). 
If we were to extend the "particle in a potential" analogy further, Eq. (E7|) could be viewed 
as the nonconservative motion of a particle. Since the rhs approaches zero in the limit 
r —>■ oo, Eq. ( p7|) gives precisely the ID asymptotic value (Eq. (p4j) for the vector potential 
A). The behavior at the origin (r = 0) is totally different; the regular singular point at the 
origin r = forces the acceptable A to vanish for m > 1 as r m . The numerical solutions of 



the two-dimensional nonlinear Schrodinger equation for m = 1, 2, 3 are shown in Fig. |TI. 
As expected, the soliton-like solutions evidently go to zero as r m for small r, and reach an 
m-independent asymptotic value predicted by Eq. ([24]) . 

More general aspects of the dynamics of the EM field can be studied mainly through 
numerical simulations of Eq. ( |PTD — this is beyond the intended scope of this paper. We 
would content ourselves here by making a few qualitative remarks and pointing out directions 
for the future efforts. The nonlinearity in Eq. (|TT| ) is a faster growing function (faster than the 
cubic) of the field amplitude, and does not exhibit saturation. We would predict then, that 
the development of the vortex chain structures from the dark stripe soliton instability for the 
general system will be faster in comparison with the low amplitude (cubic nonlinearity) case. 
Stability of vortex soliton solution is another issue that has to be dealt with. It is usually 
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FIG. 2: The form of the soliton-like solutions for the normalized potential A is shown for m = 1, 
2, 3 in case of A = 1. 

believed, though, that vortices with m = ±1 are topologically stable whereas vortices with 
larger value of the "charge" m may decay into "single-charge" vortices. 

In three dimensions (3D) the vortices form the so-called vortex line (it looks like a pancake 
with a hole in the center hanging and moving along the wire- vortex line). Effects related to 
the finite group velocity dispersion may lead to a transverse instability of the vortex line. 
All these interesting effects are left for future studies. 

We would like to emphasize that the electromagnetic fields associated with the dark 
and vortex solitons are asymptotically (at infinity) nonvanishing. Due to the generally 
accepted requirement that in physical system, the fields be localized in all directions, these 
objects have received much less attention than their localized cousins. However in recent 
experiments studying the laser field dynamics in different kinds of optical media, it has 
been demonstrated that dark and vortex solitons can be readily created as superimpositions 
upon a localized field background |J. This background can be just a few times wider than 
the soliton width. During the propagation the background spreads out, reducing its own 
intensity. In the light of these experiments let us try to put in perspective the current study 
of the dark and vortex solitons in e-p plasmas in the early Universe. Because the typical 




m=1 
m=2 
m=3 
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scale length of these solitons is the collisionless skin depth, we would need a supporting 
background spanning several skin depths. This should pose no problem because the ambient 
uniform field background could easily foot the bill. The next scale length on which we 
encounter "bullet-like" electromagnetic structures (which owe their origin to the baryon 
asymmetry) is considerably larger than the skin depth. Thus the dark and vortex solitons can 
propagate in a slowly changing background (spreading and decreasing in intensity with the 
diffractive spreading rate of the soliton decreasing as the background expands), adiabatically 
maintaining their properties, till they hit the baryon-asymmetry scalelengths. 

Topological considerations will insure the preservation of the singular points during prop- 
agation. In the propagating vortex chains, the vortices can move away from one another 
reducing the possibility of their mutual annihilation. The propagation introduces elements 
similar to Hubble expansion — the structures run away from one another. These highly 
speculative remarks need careful investigation. It is possible that the spreading of the back- 
ground field may just affect the Vortex distribution and only the cosmological expansion 
will drive them apart. 

What is extremely significant is that during the evolution of the fields, the integrals 
of motion should be preserved. It is easy to prove, by direct calculations, that Eq. ( |T7|) 
conserves the angular momentum M 



Equation (p8|) for the angular momentum is the paraxial approximation for the orbital 
angular momentum, Me — J dr[r x (E x B)], of the EM field [0. The angular momentum 
carried by the vortices is M z = mN, where N is another conserved quantity known as the 
"photon number" N = J dr ± \A ± \ 2 Q. 

It follows, then, that the relativistic e-p plasma are capable of sustaining electromagnetic 
vortex like structures, and that these structures have domains in which the EM fields carry 
non-zero angular momenta, although the total angular momentum of the entire system is 
zero. If this angular momentum could, somehow, be locally transferred to the surrounding 
medium, we would have a rather effective mechanism of imparting angular momentum to 
different domains of matter in the early universe. In our next publication we show that when 
the baryon asymmetry effects are incorporated, the medium can, indeed, acquire angular 
momentum from the EM field vortices. 




(28) 
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IV. CONCLUSIONS 



We have investigated the dynamics of the highly relativistic (7 9 3> 1) nonlinear propaga- 
tion of electromagnetic waves in unmagnetized hot electron-positron plasmas. The system 
is described by a nonlinear Schrodinger equation ( |T?D with an inverse square root type (non- 
saturating) nonlinearity. We have shown the possibility of dark and vortex soliton type 
solutions for this equation. Transverse instability of dark soliton stripes leads to the forma- 
tion of a vortex chain such that the EM fields in each vortex carry angular momentum. Such 
objects could play an important role in cosmology as sources for the structure formation in 
the MeV epoch of the evolution of the Universe. In the commonly adopted cosmological 
scenarios about the origin of the rotation of galaxies, structures grow in a hierarchy by the 
gravitational assembly of clumps out of subclumps. The origin of the angular momentum 
of galaxies if they were formed from initial fluctuations in a Friedman Universe was sug- 
gested (by Hoyle fi~8fl ) to be due to the tidal interactions between the condensing system 
Q. However, it is still not clear whether this mechanism gives an adequate solution [19|. 
We hope that the suggested mechanism of angular momentum generation in the MeV epoch 
of the Universe is an interesting alternative to explore and examine. Electromagnetism, 
operating through the versatile substrate of the e-p plasma, seems to readily generate these 
highly interesting, long-lived objects — the carriers of large amounts of mass, energy and 
angular momentum. Since an initial localization of mass, energy and angular momentum 
is precisely the seed that gravity needs for eventual structure-formation, electromagnetism 
may have provided a key element in the construction of the large-scale map of the observable 
universe. 

Results of this paper can be also applied to astrophysical objects like the pulsars, and 
active galactic nuclei — the e-p pairs are thought to be a major constituent of the plasma em- 
anating both from the pulsars, and from the inner region of the accretion disks surrounding 
the central black holes. 
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